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ABSTRACT 

We report the analysis of a highly magnetised neutron star in the Large Magellanic 
Cloud (LMC). The high mass X-ray binary pulsar Swift J045106. 8-694803 has been 
observed with Swift X-ray telescope (XRT) in 2008, The Rossi X-ray Timing Explorer 
(RXTE) in 2011 and the X-ray Multi-Mirror Mission - Newton (XMM- Newton) in 
2012. The change in spin period over these four years indicates a spin-up rate of 
—5.01 ± 0.06 s/yr, amongst the highest observed for an accreting pulsar. This spin-up 
rate can be accounted for using Ghosh and Lamb's (1979) accretion theory assuming 
it has a magnetic field of (1.2±q-7) x 10^^ Gauss. This is over the quantum critical 
field value. There are very few accreting pulsars with such high surface magnetic fields 
and this is the first of which to be discovered in the LMC. The large spin- up rate is 
consistent with Swift Burst Alert Telescope (BAT) observations which show that Swift 
J045106. 8-694803 has had a consistently high X-ray luminosity for at least five years. 
Optical spectra have been used to classify the optical counterpart of Swift J045106.8- 
694803 as a BO-1 IILV star and a possible orbital period of 21.631±0.005 days has 
been found from MACHO optical photometry. 

Key words: X-rays: binaries - stars: neutron - stars: binaries - stars: magnetar - 
individual: Swift J045106. 8-694803 



1 INTRODUCTION 

High mass X-ray binaries (HMXB) are binary systems com- 
posed of either a neutron star, white dwarf or black hole and 
an optical companion, either a supergiant star - in the case 
of supergiant X-ray binaries (SGXB) - or a dwarf, subgiant 
or giant OBe star - in the case of Be/X-ray binaries (BeXB). 
HMXB can be detected by the X-ray, optical and infrared 
emission they produce. 

X-rays are produced as matter is transferred from the 
optical star to its denser companion. It is accelerated as 
it moves into the gravitational potential well produced by 
the denser star and - if the denser star is a neutron star 
- it is then channelled by the star's magnetic field lines to 
the magnetic poles. The accreted matter then rapidly de- 
celerates when it reaches the surface and potential energy 
is converted to heat which energises the plasma, producing 
X-ray hot spots. These can appear as 'beams' if they peri- 
odically travel past our line of sight as the star rotates. The 
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pulsation period is equated with the rotation period of the 
neutron star's crust and this value changes with time as ac- 
creted matter transfers angular momentum to or from the 
star. 

HMXB are composed from massive, and therefore rel- 
atively short lived, stars and so their presence indicates 
that a new population of stars form ed relatively recently 
dCrimm. Gilfanov and SunvaevI |2003| ). This may be why 
more HMXB have been discovered in the Milky Way and 
the Small Magellanic Cloud (SMC) than the Large Magel- 
lanic Cloud (LMC). 

The HMXB in th e Milky Way are most ly found 
within the spiral arms (|Coleiro and Chatvl I2OI1I ) and the 
SMC underwent a burst of star formation about 200 mil- 
lion years ago when it was about three times closer to 
the LMC than it i s today (jCardiner and Noguchil Il996l : 
iMcBride et~ali I2OI0I ). The SMC has an HMXB popula- 
tion which is comparable in number to the Galaxy de- 
spite being almost two hundred times less massive. Whilst 
about half of the HMXB in the Galaxy are BeXB, all 
but one HMXB in the SMC is a BeXB system. The 
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Observation ID 


Start Date 


Start Time 


Exposure (ks) 


00038029001 


2008-10-23 


03:56:38 


6.56 


00038029002 


2008-11-11 


14:35:46 


6.40 


00038029003 


2008-11-14 


00:18:15 


2.59 


00038029004 


2009-01-10 


13:23:59 


4.32 


00038029005 


2011-09-30 


16:50:35 


4.47 


00038029006 


2011-10-03 


18:28:46 


5.90 



Table 1. Summary of XRT datasets. 



LMC, on the other hand, cont ains relatively few HMXB , 
most of which are BeXBs llNegueruela and Cod |2002|: 
Liu, van Paradiis and van den Hauvell l2005l : I Sturm et m 
20121 ). 

Swift J045106.8-694803 is part of a newly discovered 
BeXB system located in the LMC. It was first discov- 
ered by Beardmore et al (2009) at an RA, Dec (J2000) 
of 04:51:06.8 and -69:48:03.2 respectively, with an uncer- 
tainty of 3.5". Beardmore et al measured an X-ray flux of 
(1.68±0.11)xl0"" erg cm'^s"^ over 0.3-10 keV, flt with a 
power law of photon index 0.96±o;of with a column density 
of (1.9±0.3)xl0^^ cm"^ 

The 14-195 keV flux was found to be (2.8±0.3) x 10"" 
erg cm'^s"^ fit with a power law of photon index 2.5±0.4. 
Beardmore et al also report a pulsation period of 187 sec- 
onds. The optical companion to this X-ray source is the V = 
14.70 blue star [M2002] 9775 (Masscy 2002) with an orbital 
period of 21.64±0.02 days. 

Grebenev et al have recently calculated a photon index 
of 0.5±0.5 using data from INTECRAL. They have also cre- 
ated energy spectra over 3-200 keV showing that the energy 
of the high energy cut-off in the spectrum is at 16.0±5.0 keV 
l|Grebenev. Lutovinov and Tsvgankovll2012l '). 

In this paper we look at Swift J045106. 8-694803 in more 
detail, confirming its luminosity, pulse period, orbital pe- 
riod and spin-up rate which can be used to determine the 
strength of its magnetic field. 



2 OBSERVATIONS 
2.1 X-ray Observations 

2.1.1 Smft/XRT 

Swift's X-ray telescope (XRT) is a CCD imaging spectrom- 
eter operating over 0.2-10 keV in photon counting mode. 
Archival data were taken from NASA's High Energy As- 
trophysics Science Archive Research Center (HEASARCj^ 
covering six time periods from 23rd October 2008 to 3rd 
October 2011, as summarised in Table 1. 

The images were extracted using the /toojf] xselect. 
Source and background spectra were then extracted from 
regions of 34"radii. The spectra were binned to have 50 
counts per bin. The ancillary response files (ARF) were 
calculated with xrtmkarf and a redistribution matrix file 
(RMF) was taken from HEASARC's calibration database 
(CALDB). The position of the source was confirmed using 
ftool xrtcentroid. 

The total count rate and error of each dataset, as 

^ http://heasarc.gsfc.nasa.gov/ 
^ http:/ /heasarc. nasa.gov/ftools/ 



Observation ID 


Start Date 


Start Time 


Exposure (ks) 


96441-01-01-00 
96441-01-01-01 


2011-10-28 
2011-10-28 


13:54:55.6 
10:51:38.6 


9.005 
0.849 



Table 2. Summary of RXTE datasets. 

well as the intrinsic H I column density (Nh) and pho- 
ton index, were calculated using ftool xspec. The spec- 
tra were described by an absorbed power law with a 
fixed Galactic foreground co lumn density of 8.4x10^" cm"^ 
(Dic kev and LockmanI [199G| ) and abundances set in accor- 
dance with Wilms, Allen and McCray (2000). Intrinsic ab- 
sorption and the abun dances of elements heavier than he- 
lium were set to 0.4 (|Borkowski. Hendrick and Reynolds! 
l2006l ). X-ray spectra were then compiled in xspec over 0.2-10 
keV. 

The 0.2-10 keV fiux of each dataset was determined 
using xspec. The lumin osity was then calcul ated using a dis- 
tance of 50.6±1.6 kpc ([Bonanos et "aZl l201ll ). 

The light curve of each dataset was extracted in xselect 
and a Lomb-Scargle normalised periodogram was produced 
for each using time-series analysis package Penoqj. 

2.1.2 Swift /BAT 

The Swift Burst Alert Telescope (BAT) data were taken 
with the Hard X-ray Transient Monitor from 16th December 
2004 to the 31st May 2010 over 14-195 keV. A 58 month 
light curve was downloaded from NASA's Swift /BAT 58- 
Month Hard X-ray SurvejQ. This contained an average of 
~15 observations a day split into 8 energy bands. 

2.1.3 RXTE 

Archival data were taken from HEASARC. These were 
recorded in two datasets on the 28th October 2011 using 
the Rossi X-ray Timing Explorer's (RXTE) Proportional 
Counter Array, over ~3-10 keV, as summarised in Table 2. 

Cleaned light curves were produced for each dataset 
and combined. A Lomb-Scargle normalised periodogram was 
then created using a frequency interval of 1x10"^ Hz. 

2.2 Optical observations 

2.2.1 Optical photometry 

Archival data were taken from the MAssive Compact Halo 
Objects (MACHO) project using the 1.27 meter telescope at 
the Mount Stromlo Observatory in Australia. This covered 
the period from 1st November 1992 to the 29th December 
1999 and contains instrumental magnitudes using red (R) 
and blue (B) filters. 

The data were filtered to remove results fiagged as er- 
roneous. The four points remaining in the R band dataset 
which were over 2.4 standard deviations from the mean were 
also removed. This left 96 data points in the red and 207 in 
the blue. A Lomb-Scargle normalised periodogram was then 
created in Period for both the R and B band datasets using 

^ http: / /www. starlink.rl.ac.uk/docs/sunl67.htx/sunl67. html 
http:/ /swift. gsfc.nasa.gov/docs/swift/results/bs58mon/ 
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a frequency interval of 1x10"^ Hz and detrending using a 
polynomial of order 3. 

2.2.2 Optical spectroscopy 

Optical spectra of Swift J045106. 8-694803 have been taken 
on three separate occasions. Red end spectra were taken 
with the 1.9m Radcliffe telescope at the South African As- 
tronomical Observatory (SAAO) on the 12th December 2009 
and 26th September 2011. The data were obtained using the 
unit spectrograph combined with a 1200 lines mm"'^ grating 
and the SITe detector at the Cassegrain focus. The resulting 
spectra have a spectral resolution of ~ 3 A. Comparison cop- 
per neon spectra were taken immediately before and after 
the observation and were used for wavelength calibration. 

The source was also observed on the 8th and 10th De- 
cember 2011 with the New Technology Telescope (NTT), 
La Silla, Chile. Grisms 14 and 20 on the ESO Faint Ob- 
ject Spectrograph (EF0SC2) were used for blue and red 
end spectroscopy respectively with a slit width of 1.5". 
Grism 14 has a grating of 600 lines mm~^ that yielded 
1 A pixel" ^ dispersion over a wavelength range of AA3095- 
5085 A. Grism 20 is one of the two new Volume-Phase Holo- 
graphic grisms recently added to EFOSC2. It has a smaller 
wavelength range, from 6047-7147 A, but a superior disper- 
sion of 0.55 A pixel"^ and 1070 lines pixeP^ Filter OG530 
was used to block second order effects. The resulting spectra 
were recorded on a Loral/Lesser, thinned, AR-coated, UV 
flooded, MPP CCD with 2048x2048 pixels, at a spectral 
resolution of ~ 10 A and ~ 6 A respectively. Wavelength 
calibration was achieved using comparison spectra of He- 
lium and Argon lamps taken through out the observing run 
with the same instrument configuration. 

The data were reduced using the standard packages 
available in the Image Reduction and Analysis Facility 
(IRAF). The resulting spectra were normalized to remove 
the continuum and and a redshift correction applied corre- 
spondin g to the recession ve locity of the LMC, taken as -280 
km s"^ (jPaturel "dll2002l ). 

3 RESULTS 

3.1 X-ray Results 

3.1.1 Smft/XRT 

Figure 1 shows the positions calculated for datasets 
00038029001 (red) and 2 (blue). The first dataset has an RA, 
Dec (J2000) of 04:51:06.4 and -69:48:02.5 with a two sigma 
error radius of 3.6". The second has an RA, Dec (J2000) of 
04:51:07.0 and -69:48:03.1 with a two sigma error radius of 
3.6". These are consistent with the positions calculated by 
Beardmore et al (2009). 

The results of data extracted using xspec are sum- 
marised in Table 3. Figure 2 shows the X-ray spectra of 
dataset 00038029001 over 0.2-10 keV. This shows that the 
counts can be divided equally into four energy ranges of 0.5- 
1.5, 1.5-3, 3-4.5 and 4.5-8 keV. The X-ray luminosities are 
plotted in Figure 3, which shows that the X-ray luminosi- 
ties measured between October 2008 and January 2009 were 
about four times higher than those measured in October and 
September 2011. 




30" 



Figure 1. Uncertainty circles of the X-ray source overlaid on an 
optical background (from ESO DSS II Blue). The red circle shows 
the position found using the Swift /XRT dataset 00038029001, the 
blue circle using dataset 00038029002 and the green using the 
XMM-Newton dataset from Bartlett et al (2012, in prep). 




Energy (keV) 

Figure 2. The upper panel shows the X-ray spectrum from 
Swift/XRT dataset 00038029001 over 0.2-10 keV, here the pho- 
ton index is 0.7±0.1 with a Njj of 3±2 counts/s, the lower panel 
shows residuals. 

Figure 4 shows the Lomb-Scargle normalised peri- 
odogram for each of the six datasets. A period is only ev- 
ident in datasets 00038029001 and 2, these are 187.0±0.3 
and 186.8±0.2 seconds respectively, and show no harmonics. 
This indicates that light curves folded at these values should 
be sinusoidal. Folded fight curves for datasets 00038029001-2 
are shown in Figure 5. 

Light curves for dataset 00038029001 in each of the four 
energy ranges determined from Figure 2 are shown in Figure 

6. A plot of pulse-fraction against energy is shown in Figure 

7. The pulse-fraction appears to increase with increasing en- 
ergy between 0.5 and 4.5 keV and then decrease, however 
the large error bars make this inconclusive. 

The side peaks in the periodograms are due to gaps in 
the observations. This was confirmed firstly, by fitting a sine 
wave with the same period to the data, which gave rise to 
identical peaks. Secondly, by splitting dataset 00038029001 
into 6 shorter datasets composed of continuous observations 
and adding the individual periodograms, as can be seen in 
Figure 8. 
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Observation ID 


Count Rate 




Photon 


r iux 


Luminosity 


Reduced 


Degrees of 






(counts/s 


(10^1 cm-2) 


Index 


(10-" ergs/ 


(10^'^ ergs/s) 


Chi-Squared 


freedom 






over 0.2-10 keV) 






s/cm^) 








a 


00038029001 


0.198±0.006 


3±2 


0.7±0.1 


2.3±0.1 


7.0±0.3 


1.4 


21 


b 


00038029002 


0.175±0.005 


1±2 


0.7±0.1 


2.1±0.1 


6.3±0.4 


1.2 


17 


c 


00038029003 


0.182±0.008 


3±7 


0.6±0.2 


2.3±0.2 


7.0±0.6 


0.6 


5 


d 


00038029004 


0.181±0.006 


1±3 


0.6±0.1 


2.1±0.1 


6.5±0.4 


1.2 


11 


e 


00038029005 


0.073±0.004 


15±14 


1.1±0.4 


0.7±0.2 


2.2±0.5 


0.8 


2 


f 


00038029006 


0.058±0.003 


1±16 


0.6±0.4 


0.8±0.2 


2.3±0.5 


2.0 


2 



Table 3. Sum mary of information extracted from the spectra of Swift/XRT datasets 00038029001-6 over 0.2-10 keV using a distance of 
50.6±1.6 kpc jBonanos et aZibOllh . 
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Figure 3. Luminosity plotted against time for the Swift/XRT 
and XMM-Newton datasets over 0.2-10 keV. XMM-Newton 
dataset from Bartlett et al (2012, in prep). 



Figure 5. Light curves from Swift/XRT datasets 00038029001 
(top panel) and 00038029002 (bottom panel) over 0.2-10 keV, 
folded at 187.0 and 186.8 seconds respectively. The phase shift is 
arbitrary. 
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Figure 4. Lomb-Scargle normalised periodograms for Swift/XRT 
datasets 00038029001-6 (labeled a-f) over 0.2-10 keV. The dotted 
lines indicate the 187 second period and the two side peaks of 181 
and 193 seconds mentioned by Beardmore et al (2009). 




2.0 



Figure 6. Light curves from Swift/XRT dataset 00038029001, 
folded at 187.0 seconds, over 0.5-1.5, 1.5-3, 3-4.5 and 4.5-8 keV 
(labelled a-d). 



Datasets 00038029001 and 2 provided the best signal 
and so their hght curves were combined to determine the 
best estimation for the pulse period in late 2008. A Lomb- 
Scargle normalised periodogram yielded a maximum peak 
at 186.85±0.04 seconds as can be seen in Figure 9. Monte 
Carlo simulations give this period a 99.9% confidence level. 



Bootstrapping was conducted in order to confirm that this 
is the correct peak, and not a product of the window func- 
tion, and after 5000 iterations between 166 and 200 seconds, 
a period of 187.07±0.04 seconds was found 63% of the time 
and 181.41±0.04 seconds 37% of the time. We therefore con- 
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Figure 7. Pulse-fraction plotted against energy for Swift/XRT 
dataset 00038029001 folded at 187.0 seconds, over 0.5-1.5, 1.5-3, 
3-4.5 and 4.5-8 keV. 
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Figure 8. Lomb-Scargle normalised periodogram created by 
splitting dataset 00038029001 into 6 datasets composed of contin- 
uous observations and adding the periodograms. The dotted line 
indicates the 187.0 second period found for dataset 00038029001. 



eluded that 187.07±0.04 seconds was the most probable pe- 
riod of Swift J045106.8-694803 in October/November 2008. 

3.1.2 Smft/BAT and INTEGRAL/IBIS 

The top panel of Figure 10 shows the long-term light curve of 
Swift J045106. 8-694803 with total counts over 14-195 keV, 
binned at 28 days. The count rate appears to increase over 
time, peaking in about July 2007 and continuing to remain 
above zero for the next three years. 

This was confirmed by the light curve compiled from the 
INTEGRAL IBIS data which is shown in the second panel 
of Figure 10. The bottom panel shows the light curve of the 
optical counterpart taken from OGLE III and IV in the I- 
band. The OGLE data shows that the source was brightest 
in the I-band whilst it was barely detectable in the X-ray. 
There is a slight increase in brightness after MJD 54000 but 
this is not significant. 

A Lomb-Scargle normalised periodogram was created 



Figure 9. Lomb-Scargle normalised periodogram from 
Swift/XRT datasets 00038029001 and 2 over 0.2-10 keV 
showing a maximum peak at 186.85it0.04 seconds. The side 
peaks are due to gaps in the observations. 



for the BAT dataset using Period with a frequency interval 
of 1 X 10"^ Hz but no periods were detected. 



3.1.3 RXTE 

The Lomb-Scargle normalised periodogram is shown in plots 
(a) and (b) of Figure 11. As with the Swift/XRT data, the 
side peaks are due to gaps in the observations. This was 
confirmed by fitting a sine wave with the same period to the 
data, as can be seen in plot (c) of Figure 11. The harmonics 
at exactly 1/3 and 1/4 of the pulse period, seen in plot 
(a), indicate that the light curve should be non-sinusoidal 
and have multiple peaks when folded at the pulse period, as 
shown in Figure 12. 

The shape of the light curve contains information on 
the emission geometry from the regions close to the neutron 
stars magnetic poles. At high luminosities photons can es- 
cape from the sides of the accretion column, giving rise to 
a fan-beam pattern and a more complex profile than that 
found at lower luminosities where X-rays are generally emit- 
ted in a pencil-beam. 

The RXTE data show a pulsation period of 169.8±0.3 
seconds, 17.3±0.3 seconds less than the period evident in 
the Swift/XRT datasets taken in October and November 
2008. This gives a spin-up rate (P) of -5.8±0.1 s/yr, or 
(-1.84±0.03)xl0"^ s/s, and a spin-up time scale (Ts) of - 
30.8±0.5 yr using 



P 



(1) 



where P is the average pulse period in seconds. Since this is 
an average over three years, it is possible that the spin-up 
rate has been much higher at points during this period. 
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Figure 10. Long-term light curve from Swift/BAT (14-195 keV), INTEGRAL IBIS (15 keV - 10 MeV) and OGLE III and IV. 
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Figure 11. Panel (a) shows the Lomb-Scargle normalised peri- 
odogram from the combined RXTE datasets over 3-10 keV. Panel 
(b) shows a close up view of the region around the main peak com- 
pared to panel (c) which shows the result of a simulated dataset 
produced by a pure sine wave of period 169.78 seconds. 

3.1.4 XMM-Newton 

Swift J045106. 8-694803 was observed by the X-ray Multi- 
Mirror Mission - Newton (XM M-Newton) on 17th July 2012 
l|Bartlett et o/l2012. in prepj') . This confirmed the position of 
Swift J045106. 8-694803 - as is shown in Figure 1 - and gave 




0.0 0.5 -.0 1.5 2.0 

Phase 



Figure 12. Light curve folded at 169.78 seconds from RXTE 
datasets 96441-01-01-00 and 01 over 3-10 keV. 



an up-to-date period and luminosity. The position was found 
to be at an RA, Dec (J2000) of 04:51:06.7 and -69:48:04.2 
respectively, with a one sigma uncertainty of 1". The lumi- 
nosity was found to be (9.8±0.9) x 10'^'' ergs/s - as shown in 
Figure 3 - and the period 168.5±0.2 seconds. 

This is 1.3±0.4 seconds less than the period calculated 
from RXTE in 2011, giving a P of -1.8±0.5 s/yr. It is 
also 18.6±0.2 seconds less than the period calculated from 
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Figure 13. Plot showing the three periods measured by 
Swift/XRT in 2008, RXTE in 2011 and XMM-Newton in 2012 
(Bartlett et al, 2012, in prep). The straight line is a line of best 
fit showing that the period is continuing to decrease, despite the 
fact that the rate of change is also decreasing. 
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Figure 14. Light curves for the R and B band MACHO datasets. 

Swift/XRT in 2008, giving an average P of -5.01±0.06 s/yr 
and a Ts of -34.9±0.4 yr. Figure 13 shows the three periods 
measured by Swift/XRT in 2008, RXTE in 2011 and XMM- 
Newton in 2012. The luminosity and rate of change of spin 
period may be decreasing but Swift J045106. 8-694803 is still 
continuing to spin-up at a high rate. 

3.2 Optical results 

3.2.1 OpUcal photometry 

Figure 14 shows that the flux of the optical companion to 
Swift J045106. 8-694803 appears to have remained fairly con- 
sistent in the B and R bands over 7 years. The B band 
dataset show a possible orbital period of 21.631±0.005 days 
as seen in Figure 15, although any underlying non- radial pul- 
satio ns from the Be star may affect the results IjBird et aA 
l2012l ). Non-radial pulsations occur when some parts of the 
stellar surface move inwards while others move outwards at 
the same time. It has been suggested that they could help 



Period (Days) 

Figure 15. Lomb-Scargle normalised periodogram for the B 
band MACHO dataset showing a possible orbital period of 
21.631±0.005 days. 
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Figure 16. Light curve from the B band MACHO dataset folded 
at 21.631 days. 



form the Be stars circumstellar disk. If the rotational ve- 
locity of the Be star is close to the critical velocity, then 
pulsations lead to matte r be ing ejected and spun up to form 
a Keplerian disk (Reig 20111). 

This period is not evident in the R band, this is most 
Ukely due to the lack of data points rather than the conse- 
quence of a real effect. This idea was confirmed by randomly 
removing half of the B band data points and creating a new 
Lomb-Scargle normalised periodogram, which also failed to 
show any evidence of an orbital period. A Ught curve folded 
at 21.631 days was then produced for the B band dataset as 
shown in Figure 16. Colour ratio and colour magnitude dia- 
grams were also produced. These confirmed that the magni- 
tude in each colour band has remained fairly consistent over 
the observation period. 

Swift J045106. 8-694803 was also observed as part 
of the OGLE III (LMC136.6.14874) and OGLE IV 
(LMC531.05.4251) programmes. These I-band data cover, 
in total, a duration of more than a decade, and are shown 
in the bottom panel of Figure 10. From this figure it is clear 
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that the source undergoes a significant long-term modula- 
tion on periods in excess of 400 days. Formally, a Lomb- 
Scargle power spectrum gives the peak to be around 440 
days, but a visual inspection of the light curve shows that 
there are other time-scale changes occurring. It is therefore 
unlikely that this long period is related directly to the binary 
period, but rather may either indicate a general time-scale 
for fluctuations in the stellar wind, or precessional motion 
of the circumstellar disk. All these variations make it very 
difficult to search for confirmation of the 21.631 day pe- 
riod seen in the MACHO data. However, if only the better- 
sampled OGLE IV data are merged with the MACHO data 
(normalised to the approximate starting magnitude of the 
OGLE HI data), then the strength of the 21.631 day peak 
in the Lomb-Scargle power spectrum increases slightly. But 
without the prior knowledge of this period from the MA- 
CHO data such a period would not have been found in the 
OGLE data. 



the absolute magnitude of emission line stars, which will be 
dependent on the inclination and size of the disc. As such we 
classify the optical counterpart of Swift J045106. 8-694803 as 
a BO-1 III-V star. 

Figure 18 shows the three red end spectra of Swift 
J045106.8-694803 taken two years apart. The ESO (top) 
spectrum is offset by 6 flux units. The Ha equivalent width, 
considered an indicator for circumstellar disk size, is remark- 
ably similar in all three spectra increasing from -29±2 A 
for the SAAO spectrum taken in 2009 to -33±1 A and - 
34.5±0.6 A for the SAAO and ESO spectra taken in 2011. 
It is not uncommon to see large variations in the Ha equiv- 
alent width of these systems on timescales of months. This 
apparent consistency in the equivalent width (and hence the 
disk size) is almost certainly linked to the exceptionally per- 
sistent X-ray activity of the source. 



3.2.2 Optical spectroscopy 

OB stars in the Milky Way are classified us ing certain metal 
and helium line ratios jWalborn and Fi tzpatri cj 1990l) 
based on the Morgan-Keenan system (MK; IMorgan et aA 
However, this is unsuitable in lower metallicity en- 
vironments as the metal lines are either much weaker or not 
present. As such, the optical spectrum of IC R J05414-6858 
was classified using the method developed bv lLenno 
for B-type stars in th e SMC and i mplemented for the SMC, 
LMC and Galaxy bv lEvans etm (|2004 12006D . 

Figure 17 shows the unsmoothed optical spectrum of 
Swift J045106. 8-694803. The spectrum is dominated by the 
hydrogen Balmer series and neutral helium lines. There does 
appear to be evidence for the Hen A4200A line, but it is dif- 
ficult to distinguish above the noise level along with the Hen 
AA4541, 4686A lines. The Hei A4143A line is clearly stronger 
than the Hell A4200A line constraining the optical counter- 
part of Swift J045106.8-694803 to be later than type 09. 
There also appears to be evidence for the Siiv A4088, A4116A 
lines necessary for a Bl classification. This is supported by 
the relative strengths of the Sim A4553A and Mgii A448lA 
suggesting a classification of B2 or earlier. 

The luminosity class of the companion star was de- 
termined by the ratios of Siv A4088/Hei AA4026 - 4121 A, 
Siv A4116/Hei A4121 and Hen A4686/Hei A4713. The first 
two ratios strengthen with decreasing luminosity class (i.e. 
with increasing luminosity) whereas the latter ratio de- 
creases with increasing luminosity. The relative strengths 
of these lines are contradictory: The Hell A4686/Hel A4713 
and Siv A4116/Hel A4121 ratios suggest a luminosity class 
HI, although the proximity of the rotationally broadened 
H(5 A4102A to the Silv lines makes this more complex. The 
Siv A4088/Hel A4026 ratio is more consistent with a star 
of luminosity class V. The V band magnitude of this star 
is reported by severa l sources as between 14.6 and 14.7 
fe.g. (Massev l [2002 : Zaritskv et al 2004). An My of -4.3±0.1 
was calculated usin g a distance modulus of 18.52±0.07 
IIBonanos et a along with an mv of 14.65±0.05 and 

an extinction j4v=0.4±0.1 (calculated using the Galactic 
column density t owards the sour ce, 8.4 ± 1.0 x lO'^" cm~'^ 



4 DISCUSSION AND CONCLUSIONS 

We have confirmed the pulsation and orbital periods given 
by Beardmore et al (2009). Recent measurements show that 
Swift J045106. 8-694803 has spun-up at a rate of ~5 seconds 
a year for the last four years and although this rate is de- 
creasing with luminosity, it is still continuing to spin-up at 
a high rate. 

Swift J045106. 8-694803 may have experienced one of 
the highest spin-up rates of any known accreting pulsar but 
this can be accounted for using Ghosh and Lamb's accretion 
theory (1979) assuming that the stellar magnetic field is 
dipolar and that it is accreting from a Keplerian disk. 

Accretion theory shows that the magnetic field (B) of 
a disc-fed accreting pulsar is related to the rate of change 
of its spin period (P), its average spin period (P) and its 
average luminosity (L). 

In general, the higher the mass accretion rate (M), the 
higher the luminosity and the more angular momentum is 
accreted per second. Where magnetic stresses dominate mat- 
ter fiow at surface SI - as seen in Figure 19 from Ghosh and 
Lamb (1979) - this almost always causes spin-up. 

At surface S2, where viscous stresses dominate, this can 
add a spin-up or spin-down torque depending on the stars 
fastness parameter (oJb). This is the ratio between the angu- 
lar velocity of the neutron star and the Keplerian angular 
velocity of the disc. 

If the pulse period is long, plasma is able to penetrate 
into the magnetosphere and reach the neutron star surface, 
transferring angular momentum and causing a strong spin- 
up torque. If the pulse period is short, the plasma is unable 
to penetrate further and the neutron star spins down. 

The magnetic moment {fj) can be found if P, P and L 
are known by assuming a given mass and radius and seeing 
which value of ^ predicts the observed P using 



P = 5.0 X 1Q~^ n{ijjs)Si{M){PL%''f 



(2) 



and the resul ts oflGiiver and Ozel,200 9V This is consistent 
with B0.5III (|Wegneill2006l ). Less information is available for 



where ^30 is the magnetic moment in units of 10'''' Gauss 
cm"^ and L37 is the luminosity of the accreting star in units 
of 10^'^ ergs/s. For an 0< cjs <0.9 



n{ujs) = 1.4(1 - 2.86a;,)(l - lj,)"^ 



(3) 
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Figure 17. Spectrum of Swift J045106.8-694803 in the wavelength range AA3900-5000A with the NTT on 2011-12-08. The spectrum has 
been normahzed to remove the continuum and redshift corrected by -280 km s^^. Atomic transitions relevant to spectral classification 
have been marked. 
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Figure 18. ESO (top panel) and SAAO (middle and bottom 
panels) spectra of Swift J045106. 8-694803 in the wavelength range 
AA6400-6700A with the NTT on 2011-12-10. Spectra have been 
normalized to remove the continuum and shifted by -280 km s~^ 



within 5% accuracy and 



1.19P"'M, 



-3/7 6/7. 

7 M30 l^^j 

(4) 



5/7 



where Mn is the mass accretion rate in units of 10^''' grams/s 
and is assumed throughout to be 1.4. Si{M) and S2{M) 
are structure functions that depends on the mass equation 
of state and the dynamical response of the neutron star. 



■3/7, 



(5) 



where Re is the radius of the accreting star in units of 10® cm, 
assumed throughout to be 1, and lis is the moment of inertia 
in units of 10*^ grams cm^. 



S2{M) = (6) 
The value of the magnetic field can then be calculated using 

(7) 



R- 



Since n(ajs) depends only on cjs and ujs is inversely pro- 
portional to PLar^/'^, the P for a star of a given mass and 
magnetic moment depends only on (PLar'^''^)"^. Figure 20 
shows that Swift J045106. 8-694803 is located in the same 
region as GX 1-^4, SAX J2103. 5-^4545, 4U 2206-^54 and 
SXP 1062 on a plot of logio(P) against logio(PL37^'"'). 

GX 1-1-4 is part of a LMXB, it is a slow rotator 
with a period of about 121 seconds and has been ob- 
served to have spun-down by about 2.6 seconds a year 
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(jChakrabartv et a/lll997h . SAX J2103. 5+4545 is part of a 
HMXB system, it is also a slow rotator with a pulse period of 
358.166±0.0005 seconds and a luminosity of (2.0±0.5) x 10^® 
ergs/s. It has also been k nown to spin-up at a rate of 
4.69±0.09 seconds per year jSidoli et ffl/ll2005l ). This means 
that it has a maRnct ic field of (1.47±o:n 9) x lO^'^ Gauss. 
4U 2206+54 jR cig, T orreion and BlavlbOl^ j and SXP 1062 
jTurolla andPopovll2oil )" are also part of HMXB systems, 
both are known to be spinning down and both were recently 
identified as having extremely high magnetic fields. 

Highly magnetised neutron stars may be formed if a 
very magnetic star conserves its magnetic field during the 
neutron stars formation or, if the neutron star is rotating 
fast enough to start a dynamo effect during the first ~20 
seconds of its life. This converts heat and rotational en- 
ergy into magnetic energy and incre ases the magnetic field 
(jPerrario and Wickramasinghell200^ ). 

The first evidence for highly magnetised accreting neu- 
tron stars came from Pizzolato et al (2008), who showed 
that the X-ray source 1E161348-5055 may contain a neu- 
tron star with a surface magnetic field of ~10^^G which is 
part of a LMXB. This was shortly followed by Bozzo et al 
(2008), who argued that highly magnetised neutron stars 
may also exist in HMXB. 

Our results show that the P measured between 2008 
and 2011 implies a magnetic field of (1.3+0.1) x 10^* Gauss, 
whilst that measured between 2011 and 2012 implies a mag- 
netic field of (6+1) X 10^^ Gauss. The value calculated for 
the whole duration of observations is (1.2+0.1) x lO'^* Gauss. 
Taking into account the highest and lowest possible values of 
B measured between 2008 and 2011 and 2011 and 2012, this 
value becomes (1.2+q;7) x 10^* Gauss. Parameters calculated 
using equations 1-7 are shown in Table 4. 

All values are over the quantum critical value of 
4.4x10^^ Gauss, which means that the physics must be 
described with quantum field theory rather than classical 
physics. Although Swift: J045106. 8-694803 is powered by ac- 
cretion, rather than due to the decay and instabilities of its 
magnetic field, and hence cannot be referred to as a mag- 
netar yet, it has the potential to become one. The Galactic 



Orbital period 


21.631+0.005 days 


Pulse period in October 


187.07+0.04 s 


/November 2008 




Pulse period in October 2011 


169.8+0.3 s 


Pulse period in July 2012 


168.5+0.2 s 


Average Luminosity 


(3.4+0.3) x 10^*^ ergs/s 


P 


-5.01+0.06 s/yr 




(-1.59+0.02) X 10- s/s 


Ts 


-34.9+0.4 yr 


logio(P) 


0.70+0.01 


logio(PL37='''^) 


2.04+0.04 



n(ws) -3+f 
H (1.2+0-2) X 10^2 Gauss cm-3 
B (1.2+g-^)xl0i-^ Gauss 



Table 4. A summary of observed data and results derived using 
equations 1-7. for Swift J045106. 8-694803. 

HMXB system containing LSI 61303 recently showed mag- 
netar like behaviour when it underwent two bursts similar 
to those of soft gamma repeate rs (SGRs) jTorres et a/ll201^ : 
iPapitto. Torres and Reall2012l ). 

In conclusion, the spin-up rate of Swift J045106.8- 
694803 can be explained using Ghosh and Lamb's (1979) 
accretion theory provided that it has a magnetic field of 
(1.2+[S:?)xlO^'' Gauss. 
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the theoretical curves for three values of stellar magnetic moment assuming a neutron star mass of 1.3 M0. The dashed line shows 
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